Convection in Rotating Spherical 

Fluid Shells 

F. H. Busse* and R. Simitev 

Institute of Physics, University of Bayreuth, D-95440 Bayreuth, Germany 



"Mathematical Aspects of Natural Dynamos", E. Dormy, A.M. Soward 
(eds.), Grenoble Sciences and CRC Press, Boca Raton, pp. 119-198, ISBN- 
13:978-1-58488-954-0, 2007. 



1 Introduction 

Convection driven by thermal buoyancy in rotating spherical bodies of fluid 
has long been recognized as a fundamental process in the understanding of the 
properties of planets and stars. Since these objects are rotating in general and 
since their evolution is associated with the transport of heat from their interiors 
convection influenced by the Coriolis force does indeed play a dominant role in 
the dynamics of their fluid parts. In the case of the Earth it is the generation 
of the geomagnetic field by motions in the molten outer iron core which has 
stimulated much interest in the subject of convection in rotating spheres. But 
the zones and belts seen on Jupiter are a just as interesting phenomenon driven 
by convection in the deep atmosphere of the planet. Similarly, the differential 
rotation of the sun and its magnetic cycle are intimately connected with the 
solar convection zone encompassing the outer 29 percent of the Sun in terms 
of its radial extent. 

Theoretical studies of convection in rotating fluid spheres started about 
50 years ago. The attention was restricted to the linear problem of the on- 
set of convection and for simplicity axisymmetric motions were assumed. An 
account of these early efforts can be found in Chandrasekhar's famous trea- 
tise (1961) and in the papers by Bisshopp and Niiler (1965) and by Roberts 
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(1965). A little later it became evident that the preferred forms of convec- 
tion in the interesting limit of rapid rotation are not axisymmetric, but highly 
non-axisymmetric (Roberts, 1968). In this later paper, however, the incor- 
rect assumption was made that the preferred mode of convection exhibits a 
2;-component of the velocity field parallel to the axis that is symmetric with 
respect to the equatorial plane. The correct mode for the onset of convection 
was found by Busse (1970a) who approached the problem on the basis of the 
rotating cylindrical annulus model. This model takes advantage of the approx- 
imate validity of the Proudman- Taylor theorem and the analysis can thus be 
reduced from three to two spatial dimensions. We shall devote section 2 to 
a description of this model since it offers the simplest access to the spherical 
problem. The basic equations for the spherical problem are introduced in sec- 
tion 3 and the onset of columnar convection in spherical shells is discussed in 
section 4. In section 5 the onset of inertial mode convection is described which 
prevails at very low Prandtl numbers. In section 6 the properties of finite 
amplitude convection are outlined for moderately low Prandtl numbers Pr, 
while convection at higher values of Pr is considered in section 7. In section 
8 equatorially attached convection is considered which evolves from inertial 
mode convection. The problem of penetrative convection is addressed in sec- 
tion 9 and in section 10 some aspects of convection in the presence of thermal 
as well as chemical buoyancy are discussed. The chapter concludes with some 
remarks on applications among which the dynamo action of convection is of 
special importance. 



2 Convection in the Rotating Cylindrical An- 
nulus 

Convection in the fluid filled gap between two rigidly rotating coaxial cylinders 
is receiving increasing attention because it shares many linear and nonlinear 
dynamical properties with convection in rotating spherical fluid shells, at least 
as far as columnar convection is concerned. From the point of view of planetary 
applications it seems natural to have gravity pointing inward and to keep the 
inner cylinder at the higher temperature T2. But since the centrifugal force is 
used as a source of buoyancy in laboratory experiments with the higher tem- 
perature at the outer cylinder we shall use this latter configuration as shown 
in Fig. [H Since only the product of effective gravity and applied temperature 
gradient is physically relevant the two cases are equivalent. In experimen- 
tal realizations of the system ordinary gravity plays a minimal role when a 
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Figure 1: Sketch of the geometrical configuration of the rotating cylindrical an- 
nulus. 

vertical axis of rotation is used and when the rate of rotation is sufficiently 
high such that the centrifugal force exceeds gravity by at least a factor of two 
or three. An important ingredient of the geometrical configuration shown in 
Fig. [1] are the conical boundaries at top and bottom which cause a variation 
in height with distance from the axis of rotation. Without this variation in 
height steady two-dimensional convection rolls aligned with the axis will be 
realized since they obey the Proudman- Taylor condition. The Coriolis force is 
entirely balanced by the pressure gradient in this case and Rayleigh number 
for onset of convection in the small gap limit is given by the Rayleigh-Benard 
value for a non-rotating layer. Thin Ekman layers at the no-slip top and bot- 
tom boundaries exert only a minor infiuence on the dynamics of convection 
if the height L of the annulus is sufficiently large in comparison to the gap 
size. As soon as the height changes in the radial direction any flow involving 
a radial velocity component can no longer satisfy the geostrophic balance. In- 
stead a weak time dependence is required and the flow assumes the character 
of Rossby waves. These waves are well known in the meteorological context 
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Figure 2: The mechanism of propagation of a Rossby wave visuaUzed in the 
equatorial plane of the rotating annulus: Fluid columns originally resting at the 
mid-surface acquire anti-cyclonic vorticity relative to the rotating system when they 
are displaced outwards towards the shallow region. Cyclonic vorticity is acquired by 
the columns displaced inwards. The action of the columnar motion on the neighbor- 
ing fluid columns is such that an initial sinusoidal displacement propagates in the 
prograde direction. 

where the variation of the vertical component of rotation with latitude has 
the same effect as variation of height in the annulus of Fig. [H The dynamics 
of Rossby waves can be visualized most readily if the action of the vorticity 
acquired by the fluid columns displaced radially from the middle of the gap 
is considered. As indicated in Fig. [2] columns shifted inward acquire cyclonic 
vorticity because they are stretched owing to the increasing height. The op- 
posite sign of vorticity is exhibited by columns moving outward. Since their 
moments of inertia are increased they must rotate anti-cyclonically relative 
to the rotating system in order to conserve angular momentum. The action 
of the acquired motion of sinusoidally displaced columns on their neighbors 
results in the propagation of a wave as shown in Fig. [2l The phase velocity 
is in the prograde (retrograde) direction hen the height decreases (increases) 
with distance from the axis. 

In the case of convection the dynamics is modified by the presence of ther- 
mal buoyancy and the phase velocity is less than that of Rossby waves except 
in the limit of vanishing Prandtl number. The analysis of the thermal Rossby 
waves, as the propagating convection waves are called, is quite simple if the 
inclination of the cones with respect to the equatorial plane is introduced as a 
small perturbation. Using the gap width D as length scale and D'^/u as time 
scale where u is the kinematic viscosity of the fluid we may assume the velocity 
field in the form 



u = Vip^x, y,t) X k + . . . 
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where k is the unit vector in the 2;-direction parallel to the axis of rotation 
and where the small gap approximation with x as radial coordinate and y 
as azimuthal coordinate has been assumed. Only the geostrophic component 
of u has been denoted explicitly in expression (1). Deviations from geostro- 
phy are induced by the condition of vanishing normal velocity at the conical 
boundaries, 

L . , 

u^cVo ±u^^0 at ±— (2) 

where the tangent rjo of the angle x between the cones and the equatorial plane 
has been introduced as small parameter. By taking the 2;-component of the 
curl of the equation of motion and averaging it over the height of the anmihis 
we can incorporate the boundary condition (2) into an equation for —A21P 
which is the ^-component of vorticity, 

at ay ox ox oy J oy oy 

where A2 is the two-dimensional Laplacian, A2 = d"^ jdx^ + d"^ jdy^. Equation 
(3) must be considered together with the heat equation for the deviation © 
from the static temperature distribution of pure conduction, 

\ ot oy ox ox oy J oy 

where 6 is measured in multiples of (T2 — Ti)Pr / Ra. Ti and T2 are the 
temperatures prescribed at the inner and outer cylindrical boundaries, respec- 
tively. The dimensionless Coriolis parameter rj, the Rayleigh number Ra and 
the Prandtl number Pr are defined by 

vL UK K 

where Q denotes the angular velocity of rotation, 7 is the coefficient of thermal 
expansion, k is the thermal diffusivity and tq is the mean radius of the annulus. 
Assuming stress-free boundaries at the cylindrical walls, 

92 1 
V^ = ^V' = e = at x = ±-, (5) 

we obtain a completely specified mathematical formulation of the problem of 
centrifugally driven convection in the cylindrical annulus. 
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The onset of convection is described by the hnearized version of Eqs. (3) 
which can be solved by 

-0 — AsmmT{x + ^) ex.-p{iay + iuut), 

© = X2 , - ^ (6) 

with the following relationships for uj and it!a 

—770; 



(1 + Pr)(n%2 + Q;2)' 
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As expected the dependence of Rayleigh number on the wavenumber in the 
case of Rayleigh-Benard convection in a non-rotating layer is recovered in the 
limit 7/ = 0. The mode corresponding to n = 1 is preferred in this case, of 
course. This property continues to hold for finite 77. But as the limit 77 ^ 00 is 
approached, the values of Ra and uj do not depend on n in first approximation 
as can be seen from the following expressions for the critical values in the limit 
of large 77 



7 o -I 
12' 

,2„ 



Ra^ = 77|,(3 + 7r^77p^ + . . .) 

= -V2Pr-S|(l - ^ + ■ ■ ■) (8) 

where the definition rjp = rjPr ^/T/2{l + Pr)~^ has been used. Expressions (8) 
have been derived for 7i = 1, but they hold for arbitrary n when tt^ is replaced 
by (titt)^ The weak dependence on the radial coordinate of the problem has 
two important consequences: 

(i) The onset of convection is rather insensitive to the cyhndrical boundaries. 
The analysis can thus be apphed to the case of a sphere where these 
boundaries are missing. 

(ii) Modes of different radial dependence correspond to the same critical 
parameters asymptotically. Secondary bifurcations right above threshold 
become possible through couplings of these modes. 
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The latter possibility is indeed realized in the form of the mean flow instability. 
A transition to a solution of the form 

ip = A sai{ay — uot) sin 7r(x + -) + i? sin(a?/ — uot + Lp) sin 2tx{x + -) (9) 

occurs as the Rayleigh number is increased beyond the critical value unless the 
Prandtl number is rather small (Or and Busse, 1987). A characteristic property 
of solution (9) is the strong mean zonal shear which it generates through its 
Reynolds stress, AB sm7i{x + |)sin9?, where the bar indicates the 

average over the ^/-coordinate. Both signs of the shear are equally possible 
since the sign of B is arbitrary. The mean flow instability corresponds to a 
tilt of the convection columns as indicated in Fig. [31 When the columns are 
slightly tilted in the prograde sense towards the outside prograde momentum 
is carried outward and retrograde momentum is transported inwards leading 
to a differential rotation in which the outer fluid rotates faster than the inner 
one. An equilibrium is reached through viscous stresses which tend to oppose 
the differential rotation. The reverse situation occurs when the columns are 
tilted the other way as shown in the right plot of Fig. [31 The instability 
occurs because the differential rotation tends to increase the initial tilt and 
a feedback process is thus initiated. The mean flow instability of convection 
rolls is also possible in a non-rotating Rayleigh-Benard layer. But there it is 
usually preceded by three-dimensional instabilities. 

There is another way in which a differential rotation in the annulus can be 
generated. When curved cones instead of straight cones are used as indicated 
in Fig. [H solutions of the form (5) with separating x- and ^/-dependences are 
no longer possible. The term rjdyip in Eq. (3a) must now be replaced by 
ri{l+ex)dyil! where positive e refers to the convex case of Fig. [Hwhile a negative 
e corresponds to the concave cones of the right plot. The columns are tilted 
because the thermal Rossby wave has the tendency to propagate faster on the 
outside than on the inside when e is positive and vice versa. A differential 
rotation prograde on the outside and retrograde on the inside must thus be 
expected for e > as shown in the left plot of Fig. [H while the opposite 
results is obtained for e < 0. The experiment of Busse and Hood (1982) has 
demonstrated this effect. 

There are numerous other interesting features of convection in the cylin- 
drical annulus such as vacillations and relaxation oscillations which appear at 
higher Rayleigh numbers and which can be related to analogous phenomena 
of convection in rotating spheres. We refer to the papers by Or and Busse 
(1987), Schnaubelt and Busse (1992), Brummell and Hart (1993) for details. 
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Figure 3: The mean Bow instability leading to either outward (right figure) or 
inward (left figure) transport of prograde angular momentum. 




Figure 4: InEuence of curved conical boundaries. In the convex case (right figure) 
the columns tend to spiral outward in the prograde direction. They thus create a 
different rotation with higher angular velocity on the outside than on the inside. 
The opposite situation is found in the concave case (left figure). 
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The influence of non-axisymmetric modulations of the boundaries can also eas- 
ily be investigated in the annulus model as shown by Bell and Soward (1996), 
Herrmann and Busse (1997) and Wcstcrburg and Bussc (2003). Besides the 
narrow gap limit the finite gap case of the rotating cylindrical annulus system 
is also of interest. The discrete manifold of realizable wavenumbers gives rise 
to resonances and changes in the character of secondary instabilities. We refer 
to the recent papers by Pino et al. (2000, 2001) and Chen and Zhang (2002). 

3 Mathematical Formulation of the Problem 
of Convection in Rotating Spherical Shells 

The presence of the centrifugal force in rotating planets and stars usually 
causes some deviations from spherical symmetry. The surfaces of constant 
potential become spheroidal and a basic state of vanishing motion relative to a 
rotating frame of reference may not exist since surfaces of constant density do 
not coincide with surfaces of constant potential in general. As a result of this 
baroclinicity a differential rotation must be expected (see, for example, Busse, 
1982). But these effects are usually much smaller than those introduced by 
convection and it is thus a good approximation to neglect the effects of the 
centrifugal force and to assume that there exists a basic static solution with 
spherically symmetric distributions of gravity and temperature. 

For the theoretical description of thermal convection in rotating spheres 
usually the case of a gravity varying linearly with radius, g = —godr is assumed 
where r is the position vector with respect to the center of the sphere which is 
made dimensionless with the thickness d of the shell. We assume that a static 
state exists with the temperature distribution Ts — Tg — pd'^r'^ /2 + ATr]r~^{l — 
1])^'^ where r] denotes the ratio of inner to outer radius of the shell and j3 is 
proportional to a uniform density of heat sources in the sphere. In addition an 
applied temperature difference is admitted such that AT is the temperature 
difference between the boundaries in the special case P — 0. Of course, in 
geophysical and astrophysical applications only the super-adiabatic part of the 
temperature field must be identified with the temperature distribution given 
above. 

In addition to d, the time d'^/u and the temperature u'^/gajd'^ are used 
as scales for the dimensionless description of the problem where u denotes 
the kinematic viscosity of the fiuid and k its thermal diffusivity. We use the 
Boussinesq approximation in that we assume the density p to be constant 
except in the gravity term where its temperature dependence given by 7 = 
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{dg/dT)/g = const, is taken into account. The basic equations of motion and 
the heat equation for the deviation B from the static temperature distribution 
are thus given by 

dtu + u ■ Vu + rk X u = - Vvr + Qr + V'^u (10a) 
V ■ u = (10b) 
Pr{dtQ + u ■ Ve) = {Rtti + RaeW'^i'^ - ■ + V^O (10c) 

where the Rayleigh numbers Rai and Rae, the Coriohs number r and the 
Prandtl number Pr are defined by 

UK UK U K 

Since the velocity field it is solenoidal the general representation in terms of 
poloidal and toroidal components can be used, 

= V X (Vf X r) + Vw X r. 

By multiplying the (curl)^ and the curl of equation f llOal) by r we obtain two 
equations for v and w, 

[(V^ - dt)L2 + Td^]W^v + tQw - LaO = -r ■ V X [V X (u ■ Vu)] (11a) 
[(V^ - dt)L2 + ra^]w - rQt; = r ■ V X (n ■ Vw), (lib) 

where 5^ and denote the partial derivatives with respect to time t and with 
respect to the angle of a spherical system of coordinates r, 9, and where 
the operators L2 and Q are defined by 

L2 = -r'^V'^ + dr{r'^dr) (12a) 
Q = rcos^V^ - (L2 + r(9,)(cos^a, -r-^in^^e). (12b) 

Stress-free boundaries with fixed temperatures are most often assumed, 

V = dl,v = dr{w/r) = e = (13) 
at r = ri = rj/{\ — rj) and r = Tq = (1 — r^)^"*^. 

The numerical integration of equations (2) together with boundary conditions 
(4) in the general nonlinear case proceeds with the pseudo-spectral method as 
described by Tilgner and Busse (1997) which is based on an expansion of all 
dependent variables in spherical harmonics for the 6', ^-dependences, i.e. 



V = 

Lm 



J2vr{r,t)Pricose)exp{tm(f)} (14) 
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Figure 5: Centrifugally driven convection in a rotating fluid between an inner 

cooled and an outer heated spherical boundary. The motions are visualized through 
nearly neutrally buoyant thin platelets which align themselves with the shear. The 
photograph shows the onset of thermal Rossby waves in the form of convection 
columns in a thick fluid shell. 

and analogous expressions for the other variables, w and 0. Pf^ denotes the 
associated Legendre functions. For the r-dependence expansions in Chebychev 
polynomials are used. For further details see also Busse et al. (1998). For the 
computations to be reported in sections 5 and 6 a minimum of 33 collocation 
points in the radial direction and spherical harmonics up to the order 64 have 
been used. But in many cases the resolution was increased to 49 collocation 
points and spherical harmonics up to the order 96 or 128. 

4 The Onset of Convection in Rotating Spher- 
ical Shells 

The main difficulty in analyzing convection in rotating spherical shells arises 
from the fact that the role of the Coriolis force varies with the angle between 
gravity and the vector f2 of angular velocity. The geometrical configuration of 
the polar regions of the shell thus resembles that of a Rayleigh-Benard layer 
rotating about a vertical axis while in the equatorial region the model of the 
rotating cylindrical annulus can be applied. Only at low rotation rates does 
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convection set in in a global fashion and an axisymmetric mode can become 
preferred in this case (Geiger and Busse, 1981). At higher rotation rates the 
onset of convection does indeed occur in the form of the columnar modes as 
predicted on the basis of the annulus model of section 2. A visualisation of 
the convection motion in the form of thermal Rossby waves traveling in the 
prograde direction is shown in Fig. [51 

A rough idea of the dependence of the critical Rayleigh number -Rajc for 
the onset of convection on the parameters of the problem in the case Rae = 
can be gained from the applications of expressions (8) of the annulus model, 

4 

(PrT \ ^ 8 _i 2 

Y^^j (tan^^)3r^52-3 (15a) 



1 



2 1 



mc= ( p^ j ('^mtan6'„)S2~* (15b) 

1 



^c= (^ (i + p^)2pj 2 6(tan ^™/rj3, (15c) 

where refers to the mean radius of the fluid shell, = (r^ + ro)/2, and 
9m to the corresponding colatitude, 9m = arcsin (rm(l — ^7))- The azimuthal 
wavenumber of the preferred mode is denoted by rric and the corresponding 
angular velocity of the drift of the convection columns in the prograde direction 
is given by uJc/rric- In figure [6] the expressions (15a,c) are compared with 
accurate numerical values in the case Rae = which indicate that the general 
trend is well represented by expressions (15a,c). The same property holds 
for rric- In the case Rai = the agreement with expressions (15a, c) is not 
quite as good since the onset of convection is more concentrated towards the 
tangent cylinder touching the inner boundary at its equator because of the 
higher temperature gradient in that region. Since we shall continue to restrict 
the attention to the case Ra^. = 0, unless indicated otherwise, we shall drop 
the subscript i of Ra^. 

For the rigorous analysis of the onset problem in the limit of rapid rotation 
the dependence of the solution on the distance s from the axis must be con- 
sidered which has been neglected in the application of the annulus model. It 
turns out that a finite difference exists between the results of the local analysis 
and the exact global analysis as was pointed out already by Soward (1977). 
Using a WKBJ approach with the double turning point method of Soward and 
Jones (1983), Yano (1992) has analyzed the asymptotic problem in a refined 
version of the rotating cylindrical annulus model. He assumes a finite gap and 
the same dependence on s of the small inclination of the convex conical end 
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Figure 6: Critical Rayleigh number Raic (thick lines) and frequency ooc (right 
ordinate, thick dashed lines) as a function of the Prandtl number Pr in the case 
Tj = 0.4 for the Coriolis numbers r = 5 • 10^, 10'*, 1.5 • 10^ and 10^. The thick solid 
line corresponds to expression (15a) and (15c). 

surfaces as in the case of tlie sphere. The component of gravity parallel to the 
axis of rotation is still neglected. These assumptions have been dropped by 
Jones et al. (2000), who attacked the full spherical problem. They find that 
their results agree surprisingly well with those of Yano (1992). The analyti- 
cal findings have been confirmed through numerous numerical studies( Zhang 
and Busse, 1987; Zhang, 1991, 1992a; Ardes et al. 1997; Sun et a/., 1993) 
of which we show Fig. [7] as an example. This figure emphasizes the difference 
between the strongly spiralling nature of the convection columns at Prandtl 
numbers of the order unity or less and the more radially oriented columns at 
higher Prandtl numbers. This property is a result of the strong decrease of the 
frequency uj with increasing Pr. 

The interest in convection in rotating spherical shells has motivated a num- 
ber of laboratory investigations of the problem. The experiment of Hart et al. 
(1986) in which a spherically symmetric electric field acting on a dielectric in- 
sulating liquid is used to simulate gravity was carried out in space in order to 
avoid the interference from laboratory gravity. Less sophisticated experiments 
(Busse and Carrigan, 1976; Carrigan and Busse, 1983; Cardin and Olson, 1994; 
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Figure 7: Equatorial streamlines, rdv/dip = const., of convection columns in the 
cases Pr = 10^, r = 10^ and Ra = 3 • 10^ (upper half) and Pr = 1, r = 10^ and 
Ra = 2- 10^ (lower half). 

Cordero and Busse, 1992; Sumita and Olson, 2000) have used the centrifugal 
force with a cooled inner and a heated outer sphere to simulate the onset as 
well as the finite amplitude properties of convection. The main handicap of 
these experiments is the zonal flow generated as a thermal wind and the as- 
sociated meridional circulation in the basic axisymmetric state (Cordero and 
Busse, 1992). But this handicap can be minimized through the use of high 
rotation rates and the observations correspond quite well to the theoretical 
expectations as shown by the example of Fig. [51 

5 Onset of Inertial Convection at Low Prandtl 
Numbers 

Onset of instability in the form of inertial convection is well known from the 
case of a plane horizontal layer heated from below and rotating about a ver- 
tical axis. As has been discussed by Chandrasekhar (1961) convection in the 
form of modified inertial waves represents the preferred mode at the onset of 
instability for Prandtl numbers of the order 0.6 or less for sufficiently high 
values of the Coriolis number r. A similar situation is found for convection in 
rotating spherical shells where Zhang and Busse (1987) identified equatorially 
attached modes of convection as modified inertial waves. This connection has 
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10"= 10"^ 10° 10^ 10= 10^ lO" 



T * Pr 

Figure 8: The critical Rayleigh number Rac (solid lines) for the onset of the inertial 
convection as a function of rPr for various wavenumbers m. The dash-dotted lines 
correspond to the analytical expressions derived in Busse and Simitev (2004). For 
each m there are two modes corresponding to the two signs of to in expression (17b). 
The retrograde mode corresponding to the positive sign is preferred at lower values 
of rPr. But it turns up before the prograde mode which thus is preferred for larger 
values of rPr. 

motivated Zhang (1994, 1995) to develop a perturbation approach for analyti- 
cal description of the equatorially attached convection. Recently this approach 
has been extended and simphfied by Busse and Simitev (2004). Here we shall 
just present a short introduction to the subject. 

It is obvious from equations (11) even in their linearized versions that 
solutions for which the r— and ^—dependences separate are not admissible in 
general. Nevertheless for some parameter regimes simple, physically realistic 
solutions of the hnearized version of Eqs. (11) together with conditions (13) 
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can be obtained through the separation ansatz 



V = P^{cos9) exp{im(p + iu!t}f{r), 
w = P^_^_i{cos 9) exp{irrnp + iut}g{r), 

e = {cos 9) exp{imip + iujt}h{r). (16) 

Solutions of this form satisfy Eqs. (11) after the right hand sides have been 
dropped when the term proportional to ^^e expression for Qw can 

be neglected. This term vanishes exactly in the limit of small Pr and high r 
when the convection assumes the form of an inertial wave with the property 

r \ I r 



Jo/ Vr 

2^m(m + 2)(r/ro)-+^ 
^ (2m + l)(w(m2 + 3m + 2) -m)ro' ^ ' 



UJ 



m + 2 



(1 ± (1 + m(m + 2)(2m + 3)~^)^). (17b) 



Since this solution does not satisfy all the boundary conditions (13), weak 
Ekman layers must be added and finite critical Rayleigh numbers for onset 
have thus been obtained (Zhang, 1994). Results for the Rayleigh number for 
different values of the azimuthal wavenumber m are shown in figure [H Ac- 
cording to these results the mode with m = 1 is always preferred in the case 
?7 = if Prtau is sufficiently low (Busse and Simitev, 2004). As Prtau in- 
creases transition occurs to the mode propagating in the prograde direction. 
Besides the equatorially wall attached mode, convection can be described ap- 
proximately in the form (16) if r is less than 10^ (see Ardes et al., 1997) or if 
the thin shell limit — 1 is approached (Busse, 1970b, 1973). 



6 Evolution of Convection Columns at Mod- 
erate Prandtl Numbers 

In general the onset of convection in rotating fluid spheres occurs supercriti- 
cally. As long as the convection assumes the form of shape preserving travelling 
thermal Rossby waves as described by linear theory, its azimuthally averaged 
properties are time independent. In fact, as seen from a frame of reference 
drifting together with the convection columns the entire pattern is steady. A 
differential rotation is generated through the action of the Reynolds stress as 
explained in section 2. The latter is caused by the spiralling cross section of 
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Figure 9: Time periodic vacillations of convection at Ra = 2.8 • 10^ (left side) and 
Ra = 3 • 10^ (right side) for r = 10^, Pr = 1 The streamlines, rdv/dip = const, are 
shown in one quarter of the equatorial plane. The four quarters are equidistant in 
time (with At = 0.015 (At = 0.024) in the left (right) case in the clockwise sense 
such that approximately a full period is covered by the circles. 

the columns which persists as a dominant feature at moderate Prandtl num- 
bers far into the turbulent regime. The plots of the streamlines = const, 
in the equatorial plane shown in any of the quarter circles of figure [9] give a 
good impression of the spiralling nature of the columns. 

A true time dependence of convection develops in the form of vacillations 
after a subsequent bifurcation. First the transition to amplitude vacillations 
occurs in which case just the amplitude of convection varies periodically in 
time as exhibited in the left plot of figure [H At a somewhat higher Rayleigh 
number shape vacillations become noticeable which are characterized by pe- 
riodic changes in the structure of the columns as shown in the right plot of 
figure [91 The outer part of the columns is stretched out, breaks off and decays. 
The tendency towards breakup is caused by the fact that the local frequency 
of propagation varies with distance from the axis according to expression (15c) 
after 9m has been replaced by the local colatitude 9. 

The two types of vacillations also differ significantly in their frequencies of 
oscillation. This is evident from the time records of the energy densities of 
convection which have been plotted in figure [TOl This figure gives an overview 
of the evolution of time dependence in the interval 2.8 ■ 10^ < Ra < 10^. The 
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Figure 10: Time series of energy densities (thick solid lines), Ef (thin solid 
lines) and Nusselt number (dotted lines, right ordinate) are plotted for Pr = 1, 
T = 10"^, andRa = 2.8- 10^ 3•10^ 3.5- 10^ 7-10^ and 12- 10^ (from top to bottom). 
E^ and Ep have not been plotted. E^ is several orders of magnitude smaller than 
the other energies and Ep always approaches closely 0.4 • e( . 



various components of the energy densities are defined by 



= i(l V X (VB X r) p), Er = i(| Vffi X r 



El 



1 



V X (Vv X r) 



1 



Ww X r 



(18a) 
(18b) 
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Figure 11: Modulated shape vacillations of convection for Ra = 2.9 • 10^, r = 10^, 
Pr = 1. The plots show streamlines, = const., in the equatorial plane and are 
equidistant in time with At = 0.04 so that approximately a full period is covered. 



Figure 12: Energy densities EJ" (solid line), E'^ (thin dashed line, multiplied by 
the factor 100), E^ (short dash - long dash line), Ep (thick dashed line) and the 
Nusselt number Nu (dotted line, right ordinate) are plotted as function of Ra — Rac 
in the case r = 10^, Pr = 1. Ra^ = 1.9 • 10^ has been used corresponding to m = 10 

where v refers to the azimuthally averaged component of v and v is given by 

V = V — V. 

With a further increase of the Rayleigh number spatial modulations of the 
shape vacillations occur as shown in figure [TTl These modulations often cor- 
respond to a doubling of the azimuthal period but soon contributions with 
the azimuthal wavenumber m = 1 arise as shown in figure [TTl The pattern 
in this particular case is still periodic if appropriately sifted in azimuth. But 
as further modulations enter convection becomes quasiperiodic and with in- 
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creasing Ra a chaotic state is reached. We refer to the recent paper of Simitev 
and Busse (2003) which includes some movies to demonstrate the time de- 
pendence of convection. Figure [10] also demonstrates the diminishing fraction 
of the total kinetic energy that is associated with the poloidal component of 
motion which carries the convective heat transport. The differential rotation 
in particular increases much faster with Ra than the amplitude of convection 
since the Reynolds stress is proportional to the square of the latter. While 
this is already evident from the sequence of plots in figure [10] it is even more 
obvious from figure [I2] Here it can also be seen that the onset of vacillations 
and aperiodic time dependence tends to increase the heat transport as indi- 
cated by the Nusselt number in contrast to the situation in a planar convection 
layer with the same Prandtl number (Clever and Busse, 1987). In the latter 
case the mismatch is absent between the structure of the convection flow and 
the configuration of the boundary which inhibits the heat transport in rotat- 
ing spherical fluid shells. The time varying shift in the radial position of the 
convection columns thus promotes the heat transport. 

Surprisingly the spatio-temporal randomness of convection columns does 
not just increase at larger values of -Ra, but instead new coherent structures 
evolve (Grote and Busse, 2001). First there is localized convection as shown in 
Fig. [13] The differential rotation has become so strong that its shearing action 
inhibits convection in most parts of the spherical fluid shell as is evident from 
the figure. Only in a certain region of longitude is convection strong enough 
to overcome the shearing action of differential rotation. In the "quiet" zone 
the basic temperature profile recovers towards the purely conducting state 
and thus provides the buoyancy in the interior of the shell which sustains 
the localized convection as it is recirculated into the "active" zone by the 
differential rotation. 

After a further amplification of the differential rotation with increasing 
Ra the local intensification of convection no longer suffices to overcome the 
shearing action of the zonal flow. Instead of a spatial separation between 
"active" and "quiet" zones the system chooses a separation in time which 
manifests itself in the relaxation oscillations seen in the lowermost plot of 
Fig. [To] The fluctuating component of motion is still rather turbulent in the 
case of the relaxation oscillation as demonstrated in Figs. [TJ] and [15] When the 
differential rotation has decayed sufficiently in the near absence of Reynolds 
stresses generated by convection, a sudden burst of convection activity occurs 
leading to a sharp peak in the heat transport. But since the Reynolds stress 
grows just as suddenly as the kinetic energy of convection, the growth of the 
differential rotation occurs with only a slight delay. The shearing off of the 
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Figure 13: Localized convection for Ra = 7 • 10^,r = 1.5 • 10^, Pr = 0.5 The 
streamlines, rdv/dip = const, (first row) and the isotherms, Q = const, (second 
row), are shown in the equatorial plane for equidistant times (from left to right) 
with At = 0.03. 

convection columns then leads to their decay almost as quickly as they had 
set in. The relaxation oscillations occur over a wide region in the parameter 
space for high Rayleigh and Coriolis numbers and for Prandtl numbers of the 
order unity or less. Since it is mainly determined by the viscous decay of the 
differential rotation, the period of the relaxation oscillation does not vary much 
with these parameters. For the case of 77 = 0.4 which has been used for almost 
all numerical simulations a period of about 0.1 is usually found. 

7 Finite Amplitude Convection at Higher 
Prandtl Numbers 

The transitions from drifting convection columns to vacillating convection and 
modulated vacillating convection do not change much as as the Prandtl number 
tends to high values (Zhang, 1992b). But as Pr increases the influence of the 
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Figure 14: Relaxation oscillation of turbulent convection in the case Ra = 10^, 
r = 1.5 • lO"^, Pr = 0.5. Energy densities Et (solid line), Et (dotted line), Ep (short 
dashed line) and the Nusselt number (long dashed line, right ordinate) are shown as 
functions of time t. 




Figure 15: Sequence of plots starting at t = 0.12015 and equidistant in time 
(At = O.OlGj for the same case as in Fig. [131 Lines of constant il^ and mean 
temperature perturbation, Q = const, in the meridional plane, are shown in the left 
and right halves, respectively, of the first row. The second row show streamlines, 
rdv/d(p = const., in the equatorial plane. 
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Figure 16: Kinetic energy densities (upper left panel), Ep (upper right panel) , 

E"^ (lower left panel) and e( (lower right panel ) all multiplied by Pr^ as a function 
of the Prandtl number Pr in the case r = 5 • 10^ . The values of the Rayleigh number 

R = 5 - lo^ 6 • lo^ 8 • lo^ lo^ are denoted by empty circles and squares and full 

circles and squares, respectively. 

differential rotation which dominates the evolution of the convection columns 
for Prandtl numbers of the order unity and below diminishes rapidly. The 
feedback process exhibited most clearly by the mean flow instability discussed 
in section 2 ceases to operate at Prandtl numbers of the order 10. Above this 
value of Pr the properties of convection become nearly independent of Pr when 
the thermal time scale instead of the viscous one is used. This property which 
is familiar from Rayleigh-Benard convection in plane layers heated from below 
holds rather generally in convection systems. For this reason we have plotted in 
figure [T6] energy densities as defined by expressions (18), but multiplied by the 
factor Pr^ in order to demonstrate the tendency towards independence of Pr. 
It should be noted that most of values used in this figure have been obtained 
from dynamo computations. But the action of the Lorentz force is rather 
weak and hardly affects the Prandtl number dependence of convection. Only 
the lower left plot for the energy densities of the differential rotation shows 
the expected strong decay with increasing Pr. These energy densities do not 



23 



Figure 17: Convection in rotating spherical fluid shells in the cases r = 10^, 
Ra = A - 10^ Pr = 1 (left column) and t = 5 ■ 10^, Ra = 8 ■ 10^ Pr = 20 
(right column). Lines of constant mean azimuthal velocity are shown in the left 
halves of the upper circles and isotherms of Q are shown in the right halves. The 
plots of the middle row show streamlines, rdv/dip = const., in the equatorial plane. 
The lowermost plots indicate lines of constant Ur in the middle spherical surface, 
r = ri + 0.5. 
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Figure 18: Streamlines = const, in the equatorial plane for the case Pr = 
0.025, r = 10^ with Ra = 3.2 • 10^ 3.4 • 10^ 4 • 10^ (from left to right). 

decay to zero for Pr 0, however. The thermal wind relationship obtained 
from the azimuthal average of the 0-component of the curl of equation (10a) 

rk ■ = deQ (19) 

continues to require a finite field in the limit of high Pr. The right hand side 
is finite because the azimuthally averaged temperature field deviates strongly 
from spherical symmetry as long as the influence of rotation is significant. 

The typical differences between convection at Prandtl numbers of the order 
unity and higher values are exhibited in figure [T71 where typical results obtained 
for Pr = 1 and Pr = 20 are compared at similar values of Ra and r. The 
approximate validity of the thermal wind relationship (19) can be noticed from 
a comparison of the plots of and of G in the case Pr = 20. The convection 
columns retain their alignment with the axis of rotation with increasing Pr, 
but the spiralling nature of their radial orientation disappears. 

8 Finite Amplitude Inertial Convection 

For an analysis of nonlinear properties of equatorially attached convection we 
focus on the case Pr = 0.025 with r = 10^. The critical Rayleigh number 
for this case is Rac = 28300 corresponding to m = 10. As Ra is increased 
beyond the critical value other values of m from 7 to 12 can be realized, but 
m = 10 and lower values are usually preferred. An asymptotic perfectly peri- 
odic solution with m = 10 or m = 9 can be found only for Rayleigh numbers 
close to the critical value when computations are started from arbitrary initial 
conditions. On the other hand, perfect periodic patterns appear to be stable 
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Figure 19: Time series of energy densities of convection in the case Pr = 0.025, 
r = 10^, for Ra = 3.1 • 10^ 3.2 • 10^, 3.4 • 10^ (from bottom to top, left side) and 
3.8 • 10^, 4 • 10^, 8 • 10^ (from bottom to top, right side). The sohd, dotted and 
dashed hnes correspond to E"^ , Ef , Ep , respectively. E"^ is measured on the right 
ordinate. 

with respect to small disturbances over a more extended regime of supercrit- 
ical Rayleigh numbers. Distinct transitions like the transition to amplitude 
vacillations and to structure vacillations do not seem to exist for equatori- 
ally attached convection. Instead modulated patterns are typically already 
observed when Ra exceeds the critical value by 10% as can be seen in the 
plots of figure [181 These modulations are basically caused by the superposi- 
tion of several modes with neighboring values of the azimuthal wavenumber 
m which appear to propagate nearly independently. For example, the period 
of 1.68 ■ 10~^ visible in the energy densities shown in figure [11] in the case 
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Figure 20: Streamlines = const, in the equatorial plane for the case Pr = 
0.025, r = 10^ with Ra = 6 ■ 10^ 8 • 10^ 10^ (from left to right). 

Ra = 3.2 ■ 10^ is just about half the difference of 3.34 ■ 10~^ between the pe- 
riods of the modes m = 8 and m = 9 according the inertial wave dispersion 
relation (17b) for prograde modes. This property must be expected if a small 
component with m = 8 participates in the pattern shown in figure [18] which 
is dominated by the (m = 9)-component. The time series of energy densities 
shown in figure [T9] indicate that usually more than two modes contribute to the 
dynamics of the pattern with the exception of the case just discussed since the 
time dependence is not periodic as in the case when only two modes interact. 
The computations of the time series require a high spatial resolution together 
with a small time step. The time spans indicated in figure [19] are sufficient for 
reaching a statistically steady state of the fluctuating components of motion 
since these equilibrate on the fast thermal time scale of the order Pr~^. Only 
close to Rttc the adjustment process takes longer as can be seen in the case 
Ra = 3.1 ■ 10^ where a (m = 10)-pattern approaches its equilibrium state. The 
pattern corresponding to the other cases of figure [19] are shown in figures [TS], 
and [20] The differential rotation represented by relaxes on the viscous 
time scale and therefore takes a long time to reach its asymptotic regime in 
the examples shown in figure [19] But the differential rotation is quite weak 
such that it has a negligible effect on the other components of motion except 
in the case of the highest Rayleigh number of figure [19] Even smaller is the 
axisymmetric part of the poloidal component of motion which is not shown 
in the plots of figure [19] At higher Rayleigh numbers the convection eddies 
spread farther into the interior and in some cases become detached from the 
equator as can be seen in the plots of figure [20] In this way the convection 
eddies contribute to the heat transport from the inner boundary. But at the 
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same time they acquire the properties of the convection columns which are 
characteristic for convection at higher Prandtl numbers. Accordingly the dif- 
ferential rotation is steeply increased at Ra = 10^ and a tendency towards 
relaxation oscillation can be noticed in the upper right time series of figure [I9l 
For additional details on low Prandtl number spherical convection we refer to 
Simitev and Busse (2003). 
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